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Exploration of the small x kinematic region by the HERA experiments led to a revival of some models which 
existed before the advent of Quantum Chromo Dynamics (QCD) as the theory of the strong interactions. Pre- 
dictions of these models for the deep inelastic structure functions are compared with those given by QCD. Future 
experiments will concentrate on the large a;-region and we will discuss some issues which are important for the test 
of QCD. In particular we emphasize the next-to-next-to-leading (NNLO) order analysis of the structure functions 
and the determination of the strong coupling constant Os . We also make some critical remarks about the relevance 
of so called large corrections in the small and large a;-region. 




Figure 1. Kinematics of deep inelastic lepton- 
hadron scattering. 



1. Introduction 

Deep inelastic lepton-hadron scattering is given 
by the following process (see Fig. |l|) 

h{h) + H{p)^h{k2) + 'X'. (1) 

Here ^X' denotes any inclusive final hadronic 
state. The in and outgoing leptons are repre- 
sented by li and I2 respectively and the hadron 
is denoted by H. On the Born level the reaction 
proceeds via the exchange of one of the vector 
bosons V of the standard model which are given 
by 7, Z and W^. The kinematic variables needed 



in this paper are defined by 

2 r)2 / n P-1 
q = —Q < , V = , X — , 

^{p + qf = 2Mv -~Q^ + M'^, (2) 

where W represents the CM energy of the vir- 
tual boson-hadron system. The most interesting 
quantities under study are the structure functions 
denoted by Fi{x, Q^) (« — 2, 3, L) which provides 
us with information about the strong interaction 
between the quarks and gluons. Besides these 
structure functions one can also study other ob- 
servables for instance those which are related to 
the production of jets. However due to a lack of 
space we will limit ourselves to a discussion of the 
structure functions only. 

2. QCD and alternative models 

Immediately after the first experiments carried 
out at SLAG [ |l| |^ several models were proposed 
to explain the behaviour of the structure func- 
tions. The most prominent among them are 

a. Light Cone Expansions [ ^ 

b. Parton Model [ §| 

c. Vector Meson Dominance (VMD) [ || 

d. Regge Pole Model [ § 

^Similar experiments were also carried out at DESY see [ 
i 
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e. Dual Resonance Models 

One could also add the BFKL [ 0| approach 
which appeared later on the scene than the mod- 
els above. It can be considered as a merger of 
ideas put forward by the models mentioned un- 
der b and d. 

2.1. Light Cone Expansions 

Suppressing Lorentz indices the structure func- 
tion can be written as 



1 



47rM 



dSe^«-(p| [J(z),J(0)] (3) 



In the commutator above J{z) represents the 
electro-weak current. In the rest frame of the 
proton, with p ~ {M,0), the exponent becomes 
q.z ~ i'{zo — 2:3) — xMzs for 1/ » M. Using 
Fourier analysis the integral only gets contribu- 
tions in the region q.z = const.. Therefore the 
integrand is dominated by the space-time region 



zo - Z3 \< 



1 



Zo < 



Z3 < 



(4) 



X M X M 

Furthermore because of Einstein microcausality 
the commutator of the electro- weak currents van- 
ishes for < 0. From the latter and Eq. (Q) one 
infers that the integrand is dominated by the light 
cone region < provided is chosen to 

be sufficiently large. Therefore this variable is the 
relevant scale for the Operator Product Expan- 
sion (OPE) and the Parton Model which are the 
ingredients of perturbative QCD. However there 
are other arguments in the literature [ || leading 
to the claim that the distance between absorption 
and emission of the virtual vector boson, given 
by zq ^ Z3 — 1/xM, is the relevant scale. If 
1/xM > 2R, where R denotes the radius of the 
hadron, the virtual vector boson fluctuates into a 
hadronic system which interacts with the hadron 
H in Fig. |^ so that one deals with hadron-hadron 
collisions rather than vector boson-hadron scat- 
tering. When X is sufficiently small the inequality 
above is satisfied f\ and one enters the so called 

^ For the proton, where R = 5 GeV~^, this would mean 
that 2: < 0.1 



Regge region where 2Mi' ^ Q^. Note that at 
this moment this inequality can be experimen- 
tally only realized when the vector boson is repre- 
sented by the photon. The inequality 1/xM > 2R 
is the justification of small x physics which is at 
the basis of e.g. the Regge Pole Model and Vector 
Meson Dominance. 

When the light cone region dominates one can 
make the Operator Product Expansion (OPE) [ 



[J(2),J(0)] 



^5:c^-(zV)o^-(M^o), 

T N 



(5) 



where r and N denote the twist and spin of the 
operator O^''^ respectively. Both the coefficient 
function C^'"^ and the operator O^ '^ are under- 
stood to be renormalized and /i denotes the re- 
normalization scale which later on can be iden- 
tified with the factorization scale. The opera- 
tor product expansion has been proven for some 
quantum field theories in the context of perturba- 
tion theory. Assuming that it also holds in QCD 
one can express the A^th moment of the structure 
function as follows 



/ dxx^^^ F{x,Q^) 
Jo 



(6) 



In momentum space the operator matrix element 
and the coefficient function are defined by 

AW-^(/i2)^(p|O^--(;,^0)|p), (7) 

and 

C (^91^ = J d^z e''?-C^-^(zV') , (8) 

respectively. 

2.2. Parton Model 

For leading twist, i.e. r = 2, the results of 
the operator product expansion are reproduced 
by the Parton Model. The Bjorken scaling vari- 
able X, defined in Eq. (^), has a nice interpreta- 
tion if one describes lepton-proton scattering in 
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Figure 2. Parton Model for deep inelastic lepton- 
hadron scattering in an Infinite Momentum Frame. 



an Infinite Momentum Frame (IMF) [ g]. The 
latter is e.g. given by the following choice for the 
proton and vector boson momenta 



V={P 



2F 



(9) 



Writing the proton state | p) as a sum of multi 
parton states i.e. | p) = Yll=i I ^i) the momenta 
of the latter can be expressed as 



with 



n 

4=1 



2x,P 



(10) 



Using the IMF representation, the life time of the 
multi parton state becomes 



1 



Tpart 



2P 



Po 



En kl_^+ml 
i—1 Xi 



M2 



(11) 



If the interaction time of the vector boson with 
the partons is given by Tint ~ l/go = PjMv then 
the scattering becomes incoherent for T^art 2> Tint 
or 



1 



(12) 



The inequality above only holds if 7^ or Xi 7^ 
1. From these considerations one can derive the 
master formula for the structure function 



F(x,Q2)^^e2/,(x) 



(13) 



a=l 



where x is interpreted as the fraction of the 
proton's momentum carried away by the parton 
struck by the virtual vector boson V in Fig. 
Furthermore Cq and /a denote the parton charge 
and the parton density respectively where a runs 
over the number of active flavours in the initial 
state. Notice that in QCD the relation in Eq. jl^ ) 
will be modified after one has included radiative 
corrections. From the discussion above it is clear 
that at small x deep inelastic scattering becomes 
coherent. However the Parton Model does not 
give any information at which value of x this will 
happen. For this we have to know much more 
about QCD. Notice that the Infinite Momentum 
Frame picture is not the justification of the Par- 
ton Model. The latter can be also formulated in 
a Lorentz covariant way. 

The real justification for the light cone expan- 
sion and the Parton Model comes from QCD. 
This theory gives much more firm predictions 
when compared with the models mentioned be- 
fore and discussed hereafter. Let us enumerate 
some of them 

a. The dependence of Fix, Q^). 

b. Sum rules of the type dxF{x, Q"^) which 
only depend on the flavour group and the 
value of the strong coupling constant a^. 

c. Relations like the Callan-Gross relation [||] 
which states that the longitudinal structure 
function '--^ -I- 0{as). 

d. Factorization theorems leading to the uni- 
versality of parton densities like fa in Eq. 
(p^). This enables us to predict distri- 
butions and cross sections of other hard 
processes than deep inelastic lepton-hadron 
scattering. 

e. Rates and distributions of jets. 
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Here we want to emphasize that our theoretical 
understanding of integrated quantities, hke cross 
sections or the sum rules mentioned under b, is 
much better than what is known about differ- 
ential distributions or even about the structure 
functions. This is all related to our ignorance 
about the non-perturbative properties of QCD, 
which means that we have to rely on models or 
even worse on parametrizations to describe the x- 
dependence of the parton densities in Eq. ( p^ ) or 
the fragmentation properties of quarks and gluons 
into hadron jets. In particular the global analysis 
of the parton densities requires many arbitrary 
free parameters. In spite of the fact that these 
densities allow us to describe many hard processes 
their arbitrariness impairs the predictive power 
of perturbative QCD. It is for this reason that 
this theory has so many competitors in particu- 
lar in the small x region where diffraction plays a 
prominent role. Therefore better measurements 
of integrated quantities like the sum rules above 
are indispensable to show the validity of QCD al- 
though we realize that experimentally this is very 
hard to achieve. 

2.3. Vector Meson Dominance 

In the case that the process in Fig. |l| is 
dominated by the photon one can decompose 
it into a point-like and a hadron-like state i.e. 
I 7) =1 7pomt)+ I hadron). According to the 
parton model the hadronic state vanishes like 
l/Q^. However the argument for Vector Me- 
son Dominance (VMD) is akin to the one given 
for small a;- physics presented below Eq. (^. If 
Mv ^ (5^ + rri^ where m is the mass of the 
hadronic state the life time of the latter in the 
proton's rest frame (see Fig. ^) is given by [ |l^ 



1 



2v 



'^had 



> 2R. 



(14) 



Choosing for the hadronic state the p-meson 
(lightest vector meson) one can derive, using the 
proton radius R = 5 GeV~^, the inequality v > 
3 -|- 5Q^ which can be considered as the charac- 
teristic value of for the diffractive region. The 
inequality above corresponds to x < 0.1 which be- 
longs to the kinematic regime explored by HERA 



V 



2 R 
< > 



Figure 3. Photon q fluctuating into a vector meson 
V with momentum p' and mass m. 

p = (M ,0),g = (^.,0,0 , ^ly^ + Q^), 
p' = {^i^^ + Q^+m\ 0, 0, ^ly^ + Q^). 



but not quite by the old SLAC experiments. Be- 
fore we discuss the predictions of VMD let us first 
express the structure functions into the transverse 
axiQ^iv) and longitudinal ahiQ^^v) photo ab- 
sorption cross sections. 



Fl{x,Q' 



K 



K 
K 

^ 2M 



;{aT + (Jl) ■ 



(15) 



In the literature there are two versions of VMD. 
In the older version the hadronic state was only 
represented by the p-meson [|. Shortly after sin- 
gle vector meson dominance a second version was 
proposed called generalized VMD [ |ll|. In the 
latter case one has in principle an infinite num- 
ber of vector meson states which are produced 
in the process + e~ —^ ^ hadrons' . The most 
general properties of VMD are 



(TT(g^//) 



G{Q'^,iy)a^p{iy) , 
R{Q^,iy)aTiQ^i^). 



(16) 



In the formulae above a^p denotes the real 
photon-proton (g^ — 0) cross section and G is a 
model dependent function which can be inferred 
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from the process e"*" + e — > ^hadrons' . The ratio 
R in the two models is given by 



R 



R 



m In 



VMD, 



gen. VMD, (17) 



where mo is the averaged vector meson mass. 
Further £^{v) is defined by 



f^totjV P)\v=0 



=±1 



(18) 



where (Jtot{V p) denotes the total cross section for 
the reaction V + p ^ 'X' and Xy represents the 
helicity of the vector meson V. The main con- 
sequence of VMD is that at large we obtain 
i? ^ 1 and Ft ~ F2. This is in contrast to the 
prediction of the quark parton model which yields 
i? < 1 and Fl < i^2- The old SLAC data ruled 
out VMD with V = p because in this case ^ is 
large (i.e. ^ ^ 1.2) and R in Eq. ( p^ ) rises linearly 
in Q^. However generalized VMD works much 
better since the increase of R is logarithmic and 
the fitted value for ^ is much less (i.e. ^ — 0.171) 
than the one found for p-meson dominance. We 
have checked that the parametrisation in [ |l2| fits 
the data obtained for by the Hl-group [ Oj 
rather well in the region 10-* < a; < 5.10~ 
However at larger x values i.e. for x > 0.01 
the model in [ |l^ breaks down because both the 
transverse and longitudinal cross sections become 
negative. 

2.4. Regge Pole Model 

In the Regge Pole model the structure functions 
are parametrised as follows 



i^fc(x,Q') 



Qi(0)-1 



(19) 



Here (0) represents the Regge trajectory at zero 
momentum transfer since we deal with forward 
Compton scattering 7* + p — > 7* + p. Further 
i runs over all trajectories of the particles which 
can be exchanged between the photon and the 




1 = P,A2JJ' 




Figure 4. Regge pole exchange in forward Compton 
scattering. 



proton among which the most prominent are the 
pomeron P and the mesons A2 , f, /' (see Fig. 
H). Since the real photon-proton cross section 
behaves like cr-fpiv) ^ (zz/z/o)"'^"-'"^ a compari- 
son between Eq. ( pj| ) and Eqs. (|T^), ( p^ ) re- 
veals that VMD is a special case of the Regge 
Pole Model providing us with an explicit expres- 
sion for /3k,i{Q'^)- However the Regge Pole Model 
does not predict scaling unless one makes the ad- 
ditional assumption [ a 



so that one obtains 



.(o)-i 



(20) 



(21) 



In [ |Tj] one has obtained a good fit to the data 
for F2 in the range x < 0.07 and < < 
2000 GeV/c provided one allows for a de- 
pendence of 62,2- Moreover in addition to the soft 
pomeron P (ap(0) = 1.08) one also has to add 
a hard pomeron P' contribution (api{0) ~ 1.4). 
The assumption in Eq. ( po| ) has a curious implica- 
tion. This is revealed if we take the A^th moment 
of Eq. (pl|) which is equal to 



/ dxx^-^Fk{x,Q^ 
Jo 



(22) 



We will now show that this prediction is at vari- 
ance with perturbative QCD. In the latter case 
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the iVth moment is given by 
Jo 



(23) 

a=q,g 

Here /i represents the factorization as weU as the 
renormahzation scale. Since in the Regge Pole 
Model ai{0) is independent of the location of 
the singularities in the A^-planc (see Eq. (|2^)) 
does not depend on Q^. In the scaling parton 
model where C^^^ = 1 the Regge pole behaviour 
is also exhibited by the parton densities fa^\ 
However in perturbative QCD new poles will ap- 
pear in the iV-plane because of the behaviour 



AN) 



(N) , Q 
7a. 



(24) 



The anomalous dimensions 7^j^' contain poles in 
the A^-plane different from those given by the 
Regge Pole Model. For example we have 

(N) l_ l_ 1 



etc. 



(25) 



From Eq. (|2J) we infer that the appearance of 
the pole terms on the right hand side of Eq. (|23| ) 
depends on the value of which is in contra- 
diction with the assumption in the Regge Pole 
Model. Notice that the singularities at = 1 
disappear if the pole terms of the type 1/{N— 1)' 
are re-summed in all orders of perturbation the- 
ory using the BFKL equation [ ^. This follows 
from the BFKL characteristic function 
N- 1 



6a, 



as 



(26) 



which does not contain a singularity in jgg{N) at 
= 1. Unfortunately such an equation is not 
known for the other poles given by iV > 0. How- 
ever in principle there is no contradiction between 
perturbative QCD and the Regge Pole Model pro- 
vided one drops the assumption made in Eq. ( pO| ) 
leading to Eq. (pi]). Therefore the Regge poles 
show up in the angular momentum plane but not 
in the iV-plane (moment-plane). Notice that N 
is the quantum number associated with the angu- 
lar momentum in the four dimensional Euclidean 



plane. Hence we conclude that the assumption 
made in Eq. (po| ) (see [ |j) only works in the 
case of the old scaling parton model but not for 
QCD. Hence we disagree with the conclusions of 
[ |l^ that Regge poles contradict the singularity 
structure in the A^-plane predicted by perturba- 
tive QCD. Moreover we want to emphasize that 
the singularities at TV = 1 in 7^^^ (^), which 
are called the leading poles, pose no harm for the 
convergence of the perturbation series. First, as 
has been shown in a model [ |l^, the residues of 
the sub- leading singularities atiV^O,— 1,— 2--- 
are much larger that the one corresponding to the 
leading pole. Therefore the latter does not dom- 
inate the asymptotic behaviour of the coefficient 
function for z — s- 0. The second reason can be 
traced back to the structure function, appearing 
in Eq. (p3|), which can be written as a convolu- 
tion of parton densities and coefficient functions 
in Bjorken x-space. Since these densities fa{x/z) 
vanish in the region z x the effect of the lead- 
ing poles in the coefficient function, which behave 
in z-space like (In™ z)/z, will be considerably re- 
duced. A systematic study regarding the phe- 
nomenological meaning of these so called leading 
poles or small x-terms is presented in [ uM . 



3. Issues in QCD which are relevant for 
large x and 

3.1. NNLO analysis of the structure func- 
tions 

If one would like to have a test of perturba- 
tive QCD which can be performed on the same 
level of accuracy as achieved in e"*" physics it 
is necessary to get a next-to-next-to-leading order 
(NNLO) description of the structure functions. 
Notice that some quantities like i?(e+ 'X') 
or T{Z — > hadrons) are already known up to or- 
der a, so that they are even one order higher in 
accuracy than the NNLO expression of Fk{x,Q'^). 
The latter can be used to check the QCD predic- 
tion for its evolution with respect to and to 
extract the value of a, which can be compared 
with the results obtained from the LEP and SLC 
experiments. The ingredients for the NNLO anal- 
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ysis are the order corrected coefBcients [ 
and the order corrected anomalous dimensions 
(sphtting functions). However the three-loop con- 
tributions to the anomalous dimensions 7^^'' are 
only known for < 10 except for the leading 
rif part, where n/ denotes the number of flavours, 
which is computed in [ . For the computation 
of the integrals corresponding to the three-loop 
graphs one has to make a thorough study of the 
transcendental functions which show up in the 
splitting functions Pabiz) and the anomalous di- 
mensions jj^^^ ■ Typical functions which appear 
in Pab{z) are the Nielsen integrals corresponding 
to finite and infinite harmonic sums in 7^^''. A 
study of Nielsen integrals and harmonic sums is 
made in [ and [ ^ . The computation of 7^^'' 
for general N is very tedious and its feasibility 
is still under study. Therefore one should try to 
make an estimate of the three-loop anomalous di- 
mensions. For this estimate one can use all the in- 
formation on super-symmetric and conformal re- 
lations available in the literature [ . One also 
knows the leading and sub-leading pole terms of 
the type 7^^'* ~ ~ 1)'" from the solution 

of the BFKL characteristic function [ Q which is 
known up to NLO. One also gets some input from 
the behaviour at large N which is conjectured for 
the MS-scheme in [ According to the last 

reference the anomalous dimension behaves like 
7^^^'' InA^ in all orders of as for ^ 00. 
Another interesting question is whether in the 
MS-scheme the NNLO corrected structure func- 
tions are dominated by the coefficient functions 
rather than the anomalous dimensions so that the 
latter play a subordinate role. Using some of the 
ideas above one has already made a NNLO anal- 
ysis of the structure function F3{x,Q^) in [ p3| . 
One of the interesting results is that higher twist 
(r = 4, see Eq. (||)) contributions get smaller 
when NNLO corrections are included. Recently 
a similar analysis was performed for the singlet 
and non-singlet part of the structure function 
^2(2^7 Q^) in electro-production in [ |2^ . 

3.2. Measurement of as 

It is expected that NNLO corrections to the 
structure functions diminish the theoretical error 




Figure 5. Heavy flavour production V + q ^ q + 
Q + Q including virtual corrections {V — 7, W, Z). 
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on tts by about 50 % with respect to NLO. There 
is also a strong dependence on the cut Q^Mt ™" 
posed on the data set [ . A larger cut leads to 
a decrease of as(M^) by about 0.002. Another 
issue is the treatment of heavy flavour contribu- 
tions to the coupling constant. Usually this is 
done by imposing the following matching condi- 
tions 

as{nf,knf,ainn?) = Usijif + 1, A„^+i, am^) , (27) 

where m is the mass of the heavy flavour. In the 
literature one takes a = 1 which means that for 
> m? the heavy flavour behaves as a mass- 
less quark. However this does not follow from 
perturbation theory. An example is shown in [ 
|2^ . The heavy flavour contribution to the struc- 
ture function appears for the first time in order 
a1 (see Fig. ||) provided we assume that the ini- 
tial hadron state does not contain a heavy quark. 
The easiest way to study these contributions is to 
look at sum rules e.g. the Gross- Llewellyn Smith 
sum rule [ |28| which up to order is given by 



dx 



TT \ TT / 



-H 



{a2 



(28) 



Here ai, a2 represent the light quark and gluon 
contributions [ , [ [29|. The heavy quark contri- 
bution denoted by H [^7| behaves asymptotically 



as 



-ailn— -f 62 



TO 



(29) 



The asymptotic behaviour for m? 3> follows 
from the decoupling theorem so that at small 
scales the heavy quark does not contribute to 
the sum rule. The asymptotic behaviour for 
3> TO^ can be interpreted that the heavy 
flavour behaves like a light quark which is re- 
vealed by the mass singular logarithm. One can 



s{d) 




Figure 6. Charm production W -\- s{d) ^ c. 



check that the exact expression becomes equal to 
the asymptotic one for > 40 m? which implies 
that a in Eq. (27) should be chosen to be 40 in- 
stead of 1. Only for this value it makes sense to 
re-sum the large logarithm on the right-hand side 
in Eq. (^9|). This is achieved by absorbing the 
logarithm ln(Q^/TO^) into the running coupling 
constant 



1 



(30) 



After having removed the function H on the right- 
hand side of Eq. (|28|), as{nf,Q'^) has to be re- 
placed by as{nf + 1,Q^) so that the perturba- 
tion series for the sum rule is represented in an 
Uf + 1 flavour scheme. The lesson which one can 
draw from this calculation is that it is not the 
virtuality p'^ of the gluon, which is coupled to the 
heavy quark anti-quark pair in Fig. but the 
value of the external kinematic scale which 
is relevant for the heavy flavour threshold in the 
running coupling constant and the value of a in 
Eq. (|2^). Of course it might happen that rep- 
resents a physical observable like the momentum 
transfer in e.g. quark-quark scattering. Here it 
turns out that if —p^ > ni^ the flavour has to be 
treated as a light quark. 

3.3. Down and strange quark densities at 
large x 

One of the goals of future HERA experiments 
is to provide us with a better determination of the 
down d{x) and strange s{x) quark densities. To 
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start with the former the ratio d{x)/u{x) can be 
extracted from the cross sections of the charged 
current interactions from the ratio 



o-cc(e~ p) 



(31) 



It is not excluded that one can expect 
d{x)/u{x) — > const, when x ^ \ which is in con- 
trast with present parametrisations. Notice that 
until now information on d{x) is available from 
the following observables. 



a. The first quantity is given by 



{X,i 



(32) 



Until now the structure functions above are 
extracted from fixed target experiments so 
that one has to correct for nuclear binding 
effects [ 0. This means that d{x) might 
deviate from the parametrisations existing 
in the literature. 

Another observable is the lepton asymme- 
try in W-production in proton anti-proton 
colliders given by the reaction 



P + P 



w+ {W-) + 'X' 




Figure 7. Charm production W + s{d) 
including radiative corrections. 



(33) 



The strange quark density can be measured in 
charm quark production in the charged current 
process e^(e+) + p ^ i^e(t^e) + ''X' . The basic 
reactions are given by (Fig. |6|) 

s + W+~^c, s + W~^c. (34) 

It is found that the radiative corrections to the 
processes above are very small. The latter are 
represented by the one-loop corrections to reac- 
tions (^^ and the gluon bremsstrahlung (Fig. |^) 

s + c + g, s + W~ -fc + g. (35) 

However the determination of the strange quark 
density is hampered by the following back- 
grounds. The first one is given by the process 




Figure 8. Charm production in W + g c + s 
(gluon-boson fusion). 
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(36) 



Although this process is Cabibbo suppressed it is 
enhanced when it appears in + p — > I'e + ^X' 
because in this case the d becomes a valence quark 
and dy(x) 3> s{x). Another background is the 
gluon-boson fusion reaction (Fig. ||) 



9' 



w- 



(37) 



which dominates the cross section da / dx for x < 
0.1. Higher order corrections to the reaction 
above [ ^ are quite appreciable in particular for 
F'i{x,Q'^) but less for F2(x,Q^). Therefore one 
can only determine s{x) when the gluon density is 
accurately known in the region 5.10"'^ < a; < 0.1. 
To circumvent this problem the authors in [ |3^ 
have proposed to measure the D-mesons which 
emerge from the (anti)charm quarks. Instead of 
da / dx one has to compute 



d X d z 



(38) 



where [a = q, g) denote the heavy quark 
coefficient functions corresponding to reactions 
(||),(|5l) and (|^). The function Dc{z) describes 
the fragmentation of the c-quark into the D- 
meson with = z pc. It appears that for z < 0.2 
the gluon-boson fusion dominates charm quark 
production so that one has to impose a cut on z 
in order to suppress this background. In this case 
one obtains the quantity 



da 

d X 



dz '^^'^ 

Q 2 d X d z 



(39) 



leading to the result d a^^'-' ^ d a^'-' so that 
the Born reaction in Eq. (|3^ ) becomes dominant 
and the extraction of the strange quark density is 
possible. 

3.4. Factorization and renormalization 
scale dependence 

As one knows physical quantities like structure 
functions and cross sections are scheme indepen- 
dent which implies that they do not depend on 
the choice made for the factorization scale and 
the renormalization scale ^.r. However finite or- 
der perturbation theory violates this property as 



one can see as follows. Take as an example the 
A'^th moment of a structure function given by 



FW(Q2 



f(N) 

■'1 \ //2 



,(40) 



with the quark density (for the definition of 



see Eq. (26)) 



f{N) ^ 
J a 



as(M^) 



) _ ^(1) 



7f>/2/3o 



as{Ql)_ 
and the coefficient function 
l + a.(Q2)^(i)- 



^9 



.41) 



)] 


[as(Q')l 







(41) 



Notice that we have chosen fxp = fj,fi 



(42) 
fi for 



simplicity. Further Aq^^ is a scheme independent 
combination of the the anomalous dimensions and 
the coefficient c[^\ The latter is the order as con- 
tribution to the series expansion of the coefficient 
function Cg^'' and is scheme dependent. The low- 
est order terms in the perturbation series for the 
anomalous dimension and the beta-function are 
given by 7^^^ and (3o respectively. When we mul- 
tiply the coefficient function by the quark den- 
sity all scheme dependence cancels up to order 
fls. However the order term depends on the 
scheme and the scale fi because of the missing 
two-loop coefficient function and the three-loop 
anomalous dimension. If one varies the scale like 
what is usually done by Q /2 < fj, < 2 Q then the 
variation in is wholly due to the spurious or- 
der aj. term. Therefore one can raise the question 
whether this scale variation is a good estimate of 
the theoretical error due to the absence of higher 
order QCD corrections. Moreover the scale vari- 
ation performed on NLO quantities only probes 
the effect of the Born process on the NLO cor- 
rected cross section. The latter has the typical 
form 



~(i) 

'ab 



-(1) 
'ab 



2^cb 



(0) 



.aWlnf'^' 



11 



Here a'fj and o'^J,'' represent the lowest (Born) 
and the first order cross section respectively. The 
lowest order splitting function is denoted by P,'"^ . 

The quantity contains all first order contri- 
butions originating from new production mecha- 
nism's. From Eq. (43) one infers that all scale 



(43) 



dependent logarithms are multiplied by the Born 
cross section only. In the case that a^^^^ is large 
the latter gives a better indication about the the- 
oretical error due to higher order corrections than 
a simple scale variation of the NLO hadronic cross 
section. 



4. Re-summation of large corrections oc- 
curring at small and large x 

From the computations of splitting functions 
Pob(x) and coefficient functions Ck,a{x, (5^//i^) we 
infer that large corrections can appear at 



a. X 



(soft gluon exchanges) 



1 



(soft gluon bremsstrahlung) 



However the relevance of these corrections does 
not only depend on the quantities above but also 
on the behaviour of the parton densities. This 
one can see by a study of the structure functions 
given by 



E 



1 dz 

z 



!a{^-,/)Cu^a[z.^ . (44) 



Because of the convolution integral it might hap- 
pen that fa(x/ z, /i^) vanishes in the regions z — )■ 
or z 1 so that the total contribution to the 
structure functions is small in spite of the fact 
that the coefficient functions and the splitting 
functions blow up in these regions. Therefore a 
re-summation of these non-dominant corrections 
will lead to an overestimate of their effect on the 
physical quantities. 



Figure 9. Ladder graph in <f>g-theory. 



4.1. Large corrections at small x 

The large corrections occurring for z — > are 
due to soft gluon exchanges in the t-channel. 
They have the following form 



1 1-2 

In z 



l>2. 



In' 



l>0. 



(45) 



Because the singular behaviour is due to gluon ex- 
change it only appears in the singlet part of the 
coefficient functions and splitting functions. The 
re-summation of these logarithms is performed by 
the BFKL equation [ 0| . However one has to be 
careful with the interpretation of the solution of 
this equation which will be defined as C^^^^{z). 
One cannot claim that the latter represents the 
asymptotic form of the exact coefficient function 
QEXACTf^^-^ in the limit z ^ 0. This can be 



shown in $g-theory where the exact solution [ 33 



is known for the ladder diagrams (see Fig. ^). In 
this case one can also re-sum the leading terms 
given at each order [ ^ which can be consid- 
ered as the analogue of the BFKL-solution but 
now in the case of <I>|-theory. However this re- 
summation leads to an answer which differs by 
a factor of two and larger from the exact result 
taken in the limit z — s- 0. This means that also 
non-leading terms determine the asymptotic ex- 
pression for the ladder graphs. Since the non- 
leading order terms also appear in non-planar di- 
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agrams, which do not belong to the class of ladder 
graphs, one cannot even claim that the asymp- 
totic expression for the total coefficient function 
is determined by the ladder approximation only. 
Since the exact expression for the re-summation 
of a subset of graphs, analogous to the ladder ap- 
proximation, is unknown in QCD it is very hard 
to judge what the relation is between the BFKL 
solution and the perturbative QCD result in the 
small a;-region. 

4.2. Large corrections at large x 

The large corrections which appear in the limit 
z —^ 1 are due to soft gluon radiation. They have 
the following form when the Zth order quantities 
are presented in the MS-scheme for /i^ = 



ln^'-i(l 



1- z 



l> 1. 



1 



l>0. 



(46) 



Notice that the expression for the splitting func- 
tion is a conjecture made in [ The singular- 
ity structure given in the equation above is typ- 
ical for the non-singlet coefficient function con- 
tributing to F2 and F3. If the coefficient function 
in the above equation is convoluted with regular 
functions the singularity structure becomes much 
milder. In this case they have the form 



1) 



ai In^' ^(1 — z) 



l>2. 



(47) 



Examples of this behaviour can be found for the 
coefficient functions contributing to the longitu- 
dinal structure function and the cross section 
for direct photon production. The study of the 
logarithms which appear in the coefficient func- 
tion of Eq. ( |47| ) has a twofold purpose. 

a. Prediction of the coefficients a/ in Eq. ( p7| ) 
and a comparison with the results obtained 
from the exact expression. In this way one 
can check the exact calculation. 

b. The investigation of the coefficients a; will 
reveal how one has to re-sum the large log- 
arithms. 




Figure 10. One-loop vertex graph contributing to 

FW(g2). 



First we would like to discuss the work done in [ 
|3^ to re-sum the logarithmic terms occurring in 
the longitudinal coefRcient function contributing 
to Fl ■ Up to I — 2 one obtains the following form 
when the scales are chosen to he iip = I^-b. = Q 



.(2) 



(1), 
L,q\ 



ci';^{zA)^ci 



Cf z, 
.9 



-ln2(l-z) + (^-2C(2)) 



xln(l-z) +CaCf (C(2)-l)ln(l-z) 



— -PoCf ln(l — z) + regular terms in z , (48) 

Here Cf and Ca represent the well-known colour 
factors. The re-summation of these logarithms is 
easier when one performs the Mellin transforma- 
tion. In the limit N 00 one obtains 



C 



W,(2) 



iz,l) = 



1 



2N 
N 



(0) 

Tk 



In' 



N 



(49) 



and 7j' is the anomalous dimension of a certain 
operator determining the jet function [ The 
quantity represents the Sudakov anomalous 
dimension and it is related to the infrared struc- 
ture of a gauge theory like QCD. It also occurs 
in the vertex correction to the quark-gluon vertex 
represented by F(Cp). Using n-dimensional reg- 
ularization with e = n — 4 and p\ = p\ = ^ one 
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Figure 11. Some two-loop vertex graphs contribut- 

ing to i^(2)(Q2) 



obtains for the one-loop correction (see Fig. |T^) 





r 








L J 



(50) 



The two-loop expression (Fig. |ri|) looks like 



F(2)(g2) 



2e2 



(51) 



From the coefficients of the leading pole terms 
one infers that the re-summed Sudakov vertex 
correction gets an exponential form [ [ 
Because of the analogy between the IniV terms 
in the coefficient function and the 1/e poles in 
F{Q^) we expect that the same will happen in 
the former case. The re-summed expression for 
Eq. (H) has been found in [ |||. 

A similar attempt to re-sum the large logarith- 
mic terms has been made for the partonic cross 
sections (coefficient functions) which appear in di- 
rect photon production given by 



Hi{p^)+H2{p2)^i{Et) + 'X' , 

2 Ej" 



(52) 



In the reaction above Hi and H2 denote the in- 
coming hadrons and Et is the transverse energy 
of the photon. For this process the Ith order finite 



partonic cross section has the form 

+di^2i-i\n''\l-z) + ---] 1>1, (53) 



where a^^J denote the Born cross sections corre- 



sponding to the processes in Fig. 1 

,(0) 



Notice that 

(0) 

99 



aqg' is dominant in p — N scattering whereas ct, 
becomes more important in p — p collisions. Like 
in the case of the longitudinal coefficient function 
it is more convenient to take the Mellin transform 
of the partonic cross sections which become 



'^ab "ab 



CL21 In 



21 



+CL21-1 In 



2/- 



No 



N 
1>1. 



(54) 



The re-summation of the InA^ terms is carried 
out in [ ^ (see also [ [ ^) and leads to 
exponentiation analogous to what we have dis- 
cussed above. The goal of this re-summation was 
to get a better description of the transverse en- 
ergy Et distribution of the direct photon. It 
turns out [ ^ that the re-summation leads to 
an improvement of the factorization scale depen- 
dence with respect to the exact NLO hadronic 
cross section. However the discrepancy between 
da^^'-' /(IEt and the data, in particular those 
coming from the E706 experiment, remains. This 
occurs in the region x < 0.57 where Et < 9 GeV 
and VS — 31.6 GeV. There is also a little dis- 
crepancy for the UA6 experiment in the region 
X < 0.5 where Et < 6 GeV and y/S = 24.3 GeV. 
Notice that these regions lie outside the large x- 
region which is close to 1 where soft gluons domi- 
nate. Therefore one cannot expect that soft gluon 
re-summation will cure the discrepancy found at 
smaller x. Several proposals have been made to 
describe the data at smaller transverse energy. 
One is made in [ ^ where it is assumed that the 
initial state partons have an intrinsic transverse 
momentum. Another point of view is presented 
in [ ^ claiming that the low Et data are incon- 
sistent among the various experiments. 
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Figure 12. The Born contributions to direct photon 
production 9 + 5— >g + 7, 9 + 9— *5 + 7- 



5. Conclusions 

Summarizing this talk we conclude 

1. Since HERA started to explore the small 
a;-region one has seen a resurgence of mod- 
els which were ruled out by fixed target 
deep inelastic experiments carried out at 
large x i.e. 0.01 < x < 1. Some of them 
are discussed like for instance the Regge 
Pole Model and Vector Meson Dominance. 
The BFKL approach which became popular 
with the advent of HERA data merges ideas 
taken from the Regge Pole Model as well 
as from QCD. None of these models is able 
to predict the Q^-evolution of the deep in- 
elastic structure functions contrary to QCD 
where the evolution follows from the renor- 
malization group equations. However until 
now the cc-dependence of the structure func- 
tion cannot be predicted by QCD. There- 
fore as far as the small 3;-behaviour is con- 
cerned the aforementioned models cannot 
be ruled out since they give an equally good 
description of the data. 



2. Future experiments at HERA and fixed tar- 
get experiments at FNAL will explore the 
large x-region (x > 0.01) where there is no 
alternative model to perturbative QCD. 

3. The future deep inelastic lepton-hadron ex- 
periments have to lead to precision tests of 
QCD akin to the program carried out for 
6+ colliders like LEP and SLC. In par- 
ticular one has to perform, 

a. An accurate determination of the 
strong coupling constant which can 
be obtained from structure functions 
as well as jet distributions. 

b. A study of heavy flavour thresholds 
in as and the structure functions 

c. A better estimate of large corrections 
in the perturbation series for physical 

quantities. 

d. The computation of the three-loop 
splitting functions (anomalous dimen- 
sions) in order to give a full NNLO 
analysis of Fk{x, Q^). 
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